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Pages l6, 17, and l8: The equations beginning with equation ( 31 ) on 

page l6 and ending with equation (3^-) on page l8 should be corrected 
as follows: 
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where the equation of the lifting line is 
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When integrated by parts, equation ( 31 ) becomes 
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The potential due to a bent lifting line as shown in figure 5( a ) is 
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where m is the slope of the bent lifting line when y is positive 
(see fig. 5) • When the lifting line is used to approximate a wing, 
r( -b/2) and r(b/2) are zero. The components of the velocity for a bent 
lifting line can be found by differentiation of equation (33) . 

When dr(y 1 )/dy 1 is zero equation (33) becomes 
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By John C . Martin 

SUMMARY 


Exact and approximate methods based upon linearized supersonic flow 
theory have been developed for the calculation of the velocity potential 
and the down wash from thin wings of arbitrary plan form. Particular 
attention is given to the evaluation of the downwash in the plane of the 
wing. The applicability of the method inherently depends upon a knowledge 
of the load distribution over the plan form of the wing. General expres- 
sions for the velocity potential and downwash have been derived. Simple 
modifications of these expressions produced formulas for the velocity 
potential and downwash from arbitrary curved lifting lines. 

A complete development of all formulas, starting- with the basic 
solution of the linearized potential equation for supersonic flow, is 
given. Although the paper contains many new results, some of the results 
presented have been obtained by other methods and are given here solely 
for completeness. The general formulation of the downwash equations can 
easily be used in finding exact and approximate expressions for the other 
velocity components. 

The results of the theoretical development are used to determine the 
downwash from a pitching rectangular wing and to determine the expression 
for the local angle of attack necessary to give a specified load distri- 
bution. Comparisons of the exact and approximate values of the downwash 
for several lifting-line configurations are also presented. 


INTRODUCTION 


The calculation of the downwash fields induced by thin wings at 
supersonic speeds is necessary in order to evaluate accurately the aero- 
dynamic load distribution over the tail, an important consideration for 
structural and stability calculations. Available methods, including the 
present, are based on the linearized-time-independent flow considerations. 
These methods, in general, utilize the following well-known concepts: 
conical flows, potential doublets, vortices, and pressure doublets. The 
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conical-flow method is used in reference 1 to calculate the downwash 
flow fields in the plane of the wing for several plan forms including 
the rectangular wing. The potential-doublet method is presented in ref- 
erence 2 and is used to find the downwash behind a triangular wing over 
a range of Mach numbers in reference 3* The vortex method is developed 
for supersonic -flow applications in references 4 and 5 and is used in 
reference 6 to find the downwash due to thin wings approximated by lifting 
lines. The method used in reference 7 was the potential-doublet method; 
however, by integration by parts Ward obtained an expression for the 
velocity potential in space which agrees with the expression for the 
velocity potential in space determined herein. 

The present method is essentially a development of the pressure- 
doublet method for the calculation of downwash flow fields. This method 
leads easily to expressions for the velocity potential and the downwash 
from arbitrary curved lifting lines. These expressions in turn lead to 
approximate expressions for the downwash from lifting lines. These 
approximate expressions have many computational advantages. 

The method described herein has certain advantages in that most of 
the important results obtained by other methods can be obtained without 
difficulty by using the present approach. Integrations are performed 
only on the plan form; whereas other methods, excluding the conical-flow 
method, generally lead to integrations on the plan form and over the 
wake. An attempt wa$ made to present a fairly complete development of 
the pressure-doublet method and at the same time present results which 
have not been obtained by other methods. 

The accuracy of the developed approximate formulas is brought out 
in the evaluation of a rectangular pitching wing by use of exact (lin- 
earized) and approximate expressions. Comparisons are made between 
exact calculations of the downwash from certain lifting lines and 
approximate expressions. Other applications include some simple deriv- 
ations of known results plus an expression for the local angle of attack 
necessary to give a specified load distribution. 


SYMBOLS 


c 

b/2 

2 

P 

f 


chord 
semi span 

xp limits of integration 
pressure coefficient ^Pressure^/^pV^ 
denotes finite part of integral 


NACA TN 2135 


3 


f ( y i) 

i 

k 

hijhg 


M 


m 


n > n l> n 2 


equation of an arbitrary line 
variable index 
a constant 

yi limits of integration 

free-stream Mach number 
slope of line 
limits of summation 
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angular velocity about y-axis 
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free-stream velocity 
z = Z - z lf Xi = x - Xi , 
Cartesian coordinates 
increments in xj_ and yi, 
angle of attack 


Y i = y - Yi 
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a small positive number 
an auxiliary variable 
free-stream density 
local angle of attack 
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T 


T 


1 


0 

&0 



A0TE 


(0z) j> (0z) Ix > (0z) jjj 


area of plan form in forward Mach cone from 
point (x,y,z) 

area of plan form and wake in forward Mach cone 
from point (x,y,z) 

perturbation velocity potential 

increment in velocity potential at point 
(x,y,z) due to an elementary lifting area 
at point 

difference in partial derivative of 0 with 
respect to xp (0 XlU - 0 XlI ) 

value of potential difference at trailing edge 
(^TEu " $TEz) 

vertical perturbation velocity at point (x,y,z) 
due to regions I, II, and III, respectively 



line integral around area of plan form in 
forward Mach cone from point (x,y,0) 


Subscripts : 

u, Z upper and lower surfaces, respectively 

x,y,z partial derivatives with respect to x, y, 

and z, respectively 

TE trailing edge 


THEORY 


The theoretical development is divided into three parts: First, 

general formulas are derived for the potential and the upwash in space 
due to a thin lifting surface j the pressure distribution on the surface 
is assumed to be known. Second, because the formulas in the first part 
are in many cases difficult to evaluate, the expressions for the thin 
wing are used to develop formulas for the potential or upwash in space 
due to lifting lines. Third, because in some cases the expressions for 
lifting lines are difficult to evaluate, approximate expressions are also 
derived for the upwash from lifting lines. 
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Lifting Surfaces 

Potential in space .- The partial-differential equation satisfied by 
the perturbation velocity potential in supersonic flow is 


p2 <f0 _ ^0 _ &£ = o (1) 

ctoc^ dy2 dz^ 

A general solution of equation (l) is given in reference 8. (The boundary 
conditions for airfoils are given in reference 9«) For the purposes of the 
present paper this solution may be written in the form 

0 ( *,y,z) = (f ) u . ayi (2) 

where 0 U and 0^ are the values of 0 on the upper and lower sides of 
the surface z = Z]_. The finite part of the integral in equation (2) 
must be taken as indicated by the symbol f. The double integrals that 
arise will be dealt with according td the methods for finding the finite 
part of the multiple integrals given by Hadamard in reference 8. The 
area of integration r* is the area of the z = zi plane for which 

Xi < x - 0 ]/ y 2 + Z 2 

Equation (2) remains valid when the velocity potential is replaced 
by its partial derivative with respect to the free-stream direction. 

Since the pressure is directly proportional to 0 X , replacing 0 by 0 X 
in equation (2) has the effect of introducing a potential which is 
directly proportional to the acceleration potential. 

Since 


0(x,y,z) 





it follows from equation (2) that 
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where 


R'(X) = \/(X - xi ) 2 - |B 2 (y 2 + Z 2 ) 

The difference in is zero except on the plan form; therefore, 

the region of integration T is over the plan form only. 

It is proved in reference 4 that the order of integration of the 
finite part of a multiple integral can be changed. Equation (4) for 
zq = 0 can therefore be expressed in the form 


0(x,y,z) 





dX dxp dy^ 


(5) 


The integral 


f r x pfz 

R 3 


dX 


must be integrated from the aftercone emanating from point (xq,y]_,0) 
to the point (x,y,z) because a disturbance at point (x-^y^O) only 
affects points in the aftercone from the point (xi,y]_,0). The value of 

-zX 

the preceding integral is . 

(Y2 + z 2)r 

Equation (5) now becomes 


0(x,y,z) 


i_ rr A ^i 

2rt JJ r ( Y 2 + z 2) R 


( 6 ) 


In general, equation (6) is a finite integral; therefore, the finite- 
parts symbol has been dropped. The expression for the velocity potential 
given by equation (6) was obtained in reference 7 by a different approach 
as indicated in the introduction. Equation (6) may also be readily obtained 
from Volterra’s solution (reference 10) . 
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Upwash in space .- Perhaps the most frequently desired quantity is 
the upwash j f . It is shown in the appendix that equation (6) may be 

with respect to z under the integral signs without con- 
sideration of the variable limits. The result of differentiating both 
sides of equation (6) with respect to z is 


0z 



z 2 )r 2 + p 2 z 2 (Y 2 + z 2 jlx 
2i- dXi 

(y 2 + z 2 ) 2 r3 


dyi 


(7) 


An expression for the sidewash due to a discontinuous pressure sheet can 
be obtained in a manner similar to that used in finding equation (7) . 

Upwash in the wake . - In the z = 0 plane for points not on the plan 
form or in the wake, equation (7) reduces to the form 


0z 



A0 X , X 

dxi dyi 

y 2 r 


( 8 ) 


The restrictions on equation (8) can easily be understood by comparing 
it with equation (6) . Equation (6) is an expression for the potential 
in space due to a discontinuous pressure sheet. For points not on the 
plan form or in the wake, it can be seen from equations (6) and (8) that 


0z( x >y,O) = iim 7 0(x,y,z) 
z-*0 z 


The preceding expression shows that in order for the right side of the 
expression to be finite the potential must approach zero as z approaches 
zero. It is well-known that the potential is zero in the z = 0 plane 
for zero-thickness wings except on the plan form or in the wake where the 
potential has a finite discontinuity across this plane. Because 0 Z is 
in general finite and since the potential has a finite discontinuity 
across the z = 0 plane in the wake and on the plan form, equation (8) 
cannot hold on the plan form or in the wake. 

Equation (7) may be used to find 0 Z for points on the plan form 
or in the wake by performing the integrations for an arbitrary z and 
then setting z equal to zero. Generally the integrations for an arbi- 
trary z are very difficult. Both integrations of equation (7) must be 
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performed for points on the plan form; whereas, only the integration with 
respect to yp is needed for points in the wake. The reason that one or 

both integrations must be performed before taking the limit of 0 Z as z 

approaches zero is that the integrand tends to infinity along the line 
y 1 = y as z approaches zero. The line yp = y must therefore be 
removed from the area of integration. This line may be removed by 
removing a narrow strip, of width 6 , on each side of the line yi = y 
from the area of integration. (See fig. 1(a).) The integrations must be 
performed until the width e is removed from the limits of integration. 
The limit of 0 Z as z approaches zero may then be taken. The limit 
of 02 as z approaches zero cannot be taken until the limit of 0 Z 

as e approaches zero is taken unless the effect of the strip is con- 

sidered. Figure 1(b) shows the area of integration necessary to evaluate 
a point in the wake. It can be seen from this figure that the x 1 limits 
are independent of e; therefore, for points in the wake the xp inte- 
gration need not be performed before the limits are taken. 

The calculation of 0 Z for points on the plan form or in the wake 
can be considerably simplified by using a combination of equation (7) and 
equation (8). For points in the wake, divide the area of integration as 
shown in figure 2. Each region of the integration may be considered as 
an independent wing with a local angle of attack necessary to give the 
load distribution of the region. Since regions I and III can be con- 
sidered as being two different plan forms, equation (8) can be applied 
to regions I and III, but equation (T) must be applied to region II. 

The value of e is chosen small; therefore, in region II the variation 
of A0 X with yi can be neglected in the first approximation provided 
the pressure is continuous over the plan form. Pressure discontinuities 
may be dealt with by subdividing the regions of integration. 

The contribution of regions I and III to the vertical perturbation 
velocity when z = 0 is 


(0z) 


i4x i x 


dx-L dy x 


(9) 


I+III 2 *JJ I+ iii Y 2 R 
For a point above the wake the effect of region II is approximately 


<k>n s 
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where A0 x (xi,y,O) is the difference in 0 X across the lifting surface 
along the line yi = y. The xp limits are denoted by Lp and Lg- 
After integration with respect to yp equation (10) becomes 

„ i_ r L 2 A0 Xl (xi,y,O)XY(R2 - (3 2 z 2 ) 

11 2 n J Li (x 2 - p 2 z 2 )^Y 2 + z 2 )r 

and for z = O 


y+e 

dx-^ 

y-e 




^(xp^O) 


r 


- s 2 € 2 


X 


dxp 


( 11 ) 


For points in the wake the resultant vertical perturbation velocity 
is then given by the sum of equations (9) and (11): 


0 z (x,y,°) 


A0; 


XI 


2n JJ I+III Y 2 R 


dx n 


iyi - 


1_ 

Jte 


r L 2 
'll 




(xi,y,0) \Ji C 2 - P 2 ^ 2 
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( 12 ) 


Equation (12) is approximate when € is not zero; however, the accuracy 
of the approximation increases as € approaches zero and of course in 
the limit the equation is exact. The exact vertical perturbation velocity 
is then given by # 
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Equation (13) can be expressed in the form 


0 z ( x >y>°) 


i 

“Sn 


pn A 0xiYl R 

JJt ~ 


dXj_ 


ayp 



(HO 


where the second integral is a line integral around the region of the 
plan form affecting the point (x,y,0). Equation (l4) was obtained from 
equation ( 13 ) *by integrating the first term by parts with respect to yp. 

A consideration of the previous division of the plan form into 
regions with regard to vortex distributions is helpful in understanding 
the physical meaning of the preceding manipulations. 

Equation (10) is the expression for the upwash from a series of 
horseshoe vortices distributed over region II. The y^ component of the 

vortex strength is proportional to A0 x (xp,y,O) . The xp component of 

the vortex strength is zero except along y - € and y + e. (See fig. 3*) 
The spanwise distribution of the trailing vortices consists of only 
two vortices of finite strength A0(x,y,O) located along 


yp = y - € 


and along 


yi = y + e 


Equation (9) can be considered as the expression for the vertical 
perturbation velocity due to a system of vortices. The integral over 
region I represents the effect of the bound and trailing vortices 
associated with region I plus a finite vortex (see fig. 3) along 


yp = y - c 


This finite vortex is the sum of the y-components of the vortex strength 
along 


yp = y - € 
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Since the divergence of a vortex field is zero, it is not surprising 
that a vortex of finite strength exists along this line. The integral 
over region III represents the effects of the hound and trailing vortices 
associated with region III plus a vortex of finite strength along 


yi = y + € 


Equation (9) gives the effect on the field point of the vortex distri- 
bution on regions I and III, the trailing vortices associated with the 
vortices on regions I and III, and two vortices of finite strength 
located a distance e on either side of the line yi = y. For any small 
value of e the strengths of the two finite vortices are almost equal 
and, since the field point lies midway between these finite vortices, it 
follows that equation (9) should tend to infinity as e approaches zero. 
The infinity arising from equation (9) as s 0 is exactly canceled by 
the infinity arising from equation (ll) as e 0. A finite downwash 
then results at the point (x,y,0) within the wake region provided that 
the spanwise derivative of the load distribution is continuous at y. 


Upwash on the plan form .- In order to simplify the calculation of 0 Z 
on the plan form, consider an arbitrary point slightly above the plan 
form and divide the area of integration into regions as shown in figure 4. 
Region IV is small; therefore, the value of A0 Xl is approximately 

constant in region IV. For points slightly above the plan form the 
contribution of region IV to the vertical perturbation velocity is 
approximately 


M x ~A0 x (x,y,O) 
^ z 'lV St 





z^)r^ + |3^z^(y^ + 
( Y 2 + Z 2)2r3 


dx-L dyp (15) 


The result of performing the integrations in equation (15) is 


(0z)jV ~ "|- A0 x (x,y,O) (1 6 ) 

Actually the point at the apex of the hyperbola is a singular point and 
must be removed from the region of integration by a method such as is given 
in reference 8, pages 1 V 7 to 150. The integrals over regions I, II, and III 
for points slightly above the plan form are given approximately, by equa- 
tion (l4) # The vertical perturbation velocity 0 Z is then approximately 
given by (equation (l4) plus equation ( 16 ) ) : 
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2« Jmc XY 


dx]_ 


- | ^ x (x,y,°) 


(17) 


Equation (17) is independent of e; therefore, it is exact for points on 
the plan form. The result given by equation (17) was also obtained 
independently by Dr. A. Busemann of the Langley Laboratory in an unpub- 
lished analysis using a different approach. 


Lifting Lines 

Integrations that arise in evaluating the potential or the downwash 
from thin wings are in most cases very difficult to perform. For most 
downwash problems reference 6 indicates that a lifting line can be used 
as a very good approximation. In this reference, formulas are developed 
for the upwash in space due to lifting lines. In the present paper, 
formulas are developed for the velocity potential and the upwash in space 
due to lifting lines by using the same approach that was used to develop 
the previous formulas for the thin wings. The results for the upwash 
agree with the corresponding results of reference 6. 

Potential due to a lifting line .- The potential or upwash due to an 
arbitrary curved lifting line can be found from the preceding results. 

The following expression for the infinitesimal increment in the potential 
due to an elementary lifting area can be obtained from equation (6): 


50(x,y,z) _ A 0xi zX 
5xi & Yl 2it(Y 2 + z 2 )R 

The potential from an element of a lifting line is then 


( 18 ) 


cL0(x,y,z) 


A0xi 8x l zX 5 Yi 

lim 

6x l8yi 0 2ji(y2 + z2)e 


(19) 


where the product A0 X -^ Bxi is held constant as Sxi approaches zero. 

The product A0 X ^ Sx-^ is the difference in potential because A0 X -^ is 

constant in the x-direction. Since A0 is also the circulation, equa- 
tion ( 19 ) can he written 
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d0(x,y,z) = — — — dy-, (20) 

2n(Y 2 + z 2 )R 

where T denotes the circulation at point (xi,yp). 

Let the equation of the curve denoting an arbitrary lifting line be 


xi = f (yi) 


( 21 ) 


Substituting this function into equation (20) and integrating yields the 
equation 


0(x,y,z) = — 


‘’ h 2 


z(x - f)r(yi) 


*1 (Y 2 + z 2 )\/(x - f) 2 - p 2 (Y 2 + z 2 ) 


dy-L (22) 


Upwash due to a lifting line .- The derivative of equation (22) with 
respect to z is 



[(x- f) 2 - p 2 Y 2 ] 

I* 2 - 

* 2 I 

[(x - f ) 2 - p 2 (y 2 + 2z 2 )]} (x - f )r 

z ’ ’ 2n 

(Y 2 + : 

> 2 ) 2 \ 

[(x 

- f) 2 - P 2 (Y 2 + z 2 )] 

p/2 


(23) 


An expression for the sidewash due to a lifting line can be obtained in 
a similar manner. For points in the z = 0 plane not directly behind 
the lifting line, equation (23) reduces to 


0 z (x,y,O) 


— f h2 (x - f)r(yi) ^ 
2 ,T ^h 1 y 2 \j(x - f ) 2 - p 2 Y 2 1 


(24) 


For points directly behind the lifting line, however, the equation of the 
lifting line must be known in order no evaluate equation (23) in the 
z = 0 plane. 
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Unbent lifting lines.- For a lifting line parallel to the y-axis 
located at x = a, equation (23) becomes 


0z(*,y,z) = 



[(x- a ) 2 - P^^Y 2 - 

- 2 ! 

[(x - a ) 2 - p 2 (Y 2 + 

OJ 

tsl 

OJ 

[} (x - a)r 

: J h]_ 

^Y 2 + z 2 

T\ 


a ) 2 - p 2 (Y 2 + z 2 )J 

3/2 


dy l 


(25) 


When integrated by parts, equation (25) becomes 


0z( x >y> z ) — 


xy(r 2 - p 2 z 2 )r(yi) 


2it(x 2 - p 2 z 2 ) (Y 2 + Z 2 )R 
h 2 


h 2 


h 


XY(R 2 - p2 z 2) dr ^ y ^ 
dyp 


hp (X 2 - p 2 z 2 ) (Y 2 + z 2 )R 


dyp 


(26) 


When one or both limits of the first term of equation (26) are the inter- 
section of the Mach cone from point (x,y,z) with the lifting line, these 
limits are neglected because only the finite part is to be taken. For 
lifting lines that represent airfoils, T is zero at each end of the line 
therefore, in this case equation (26) becomes 


0 z (x,y>z) = -r- 


' h 2 


XY(R 2 - p 2 z 2 )^^ 


2jt > 


h l (X 2 - p 2 z 2 ) (Y 2 + z 2 )R 


dyi 


( 27 ) 


This result was obtained in reference 6 by using vortex theory. 


dT(y x ) 

dy^ 


is zero, equation (26) becomes 


When 
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0 z (x,y,z) = 


xr 


2tf(X 2 - P 2 z 2 ) 



l x2 - e 2 1 

l*- 

■ h 2 ) 2 + 2z 

-i 

OJ 

A 

1 

& 

- h 2 ) 2 + 

* 2 ] 

\Jx 2 - (3 2 1 

[<* 

r CVJ 1 
Csl 

+ 

OJ 

1 


-|x 2 - (3 2 [(y - h x ) 2 + 2z 2 ]j(y - h x ) 

[(y - hi) 2 + z 2 ] \^ 2 - (B 2 [(y - hq) 2 + z 2 ] 


(28) 


Equation (28) is the upwash from a horseshoe vortex; this result was first 
obtained in reference 11. 

The velocity potential in space due to an unbent lifting line can be 

d. r 

expressed in terms of r and . Integrating equation (22) by parts 

dyp 

gives 


a, \ r(yq) _i zR 
0(x,y,2) - —— tan - 


h2 i f ha t -1 ZR dr , 




dr 


For a horseshoe vortex is zero and 

dyq 


(29) 


it hi) = r(h 2 ) 

therefore, the potential in space due to a horseshoe vortex is 


0(x,y,z) = f 
2« 


tan 


_1 z \ju - a) 2 - (B 2 (y - h 2 ) 2 - 3 2 z 2 


(y - h 2 )X 


tan-i z \/(x - a) 2 - p 2 (y - hq) 2 - p 2 z 2 

(y - h]_)x 


( 30 ) 


The components of the velocity can he found from equation (30) by 
differentiation. 
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Lifting line of constant slope.- The potential due to a lifting line 

of constant slope can he expressed as (from equation (22)) 



where the equation of the lifting line is 


y 1 = mx 1 + k 

When integrated by parts, equation (31) becomes 



r h 2 dr(yx) 

J h-L d ^l 


tan' 



p 2 (Y 2 + z 2 ) 
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The potential due to a bent lifting line as shown in figure 5(a) is 
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1 I dT . _i M \ m / 

I dyT tan 7 ^ 


J -b/2 


y x - ^-J_ k + z i 


m / m 


pb/2 


"J ^ tan-1 

2* dy x 

Jo 


‘V* 


yi - k y 

in / 
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(33) 
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where m is the slope of the bent lifting line when y is positive 
(see fig. 5). When the lifting line is used to approximate a wing, 
r(-t/2) and T(b/2) are zero. The components of the velocity for a bent 
lifting line can be found by differentiation of equation (33). 

When dr(y 1 )/dyp is zero equation (33) becomes 





(34) 


Eauation (34) is the expression for the potential in space due to a 
vortex of constant strength as shown in figure 5(b). Since the velocity 
components can be found by differentiation, the upwash is given by 
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(35) 


where the origin of the axes has been transformed to the point of the original 

system of axes. Equation (35) can be obtained from the results of reference 6. 
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Approximations to the Upwash from Lifting Lines 

Although the integrals that arise when lifting lines are used are 
much easier to evaluate than the integrals of exact lifting-surface 
expressions, many of the integrals that arise in connection with the 
use of lifting lines cannot he evaluated in closed form and many 
involve singularities. It is therefore desirable to obtain approximate 
expressions for lifting lines. 

Approximation of an unbent lifting line .- Equation ( 28 ) can be used 
to approximate an unbent lifting line by a series of horseshoe vortices. 
The upwash from a series of horseshoe vortices that approximate the 
circulation by a series of step ^ as shown in figure 6 (a) is 




n 



i=0 


XY t (x2 - P^Yj 2 - 2pg Z g)|r(yi) - rfy-^)] 
(X 2 - p2 z 2)( Yi 2 + z 2) \J X 2 _ ^ ± 2 _ p 2 z 2 


(36) 


where Yp denotes y - yp and i takes on all integral values from 0 
to n. The upwash from a series of horseshoe vortices that approximates 
the circulation by a series of steps as shown in figure 6 (b) is 


1 V - p 2 Y i 2 - gp 2 z 2 ) [r(yn-i) - r(yj)] 

2,1 l— (X 2 - P 2 z 2 )(Y ± 2 + z 2 ) ^X 2 - P 2 ?! 2 - (3 2 z 2 


(37) 


The average of equations ( 36 ) and (37) is 


0 z 


n 



i=0 


XYj (x 2 - pgYj 2 - 2p 2 z 2 ) [r(y 1+1 ) - r(yj-i)] 
(X 2 - (3 2 z 2 )(Yi 2 + z 2 ) \/x 2 - p 2 (Y ± 2 + z 2 ) 


(38) 


Equation ( 38 ) is the upwash from a series of horseshoe vortices which 
approximates the circulation by a series of steps as shown in figure 6 (c). 
In general, for points directly behind the lifting line, equation ( 38 ) 
should give best agreement with the exact lifting line when y (the 
coordinate for the field point) is midway between consecutive values of 
y i . An expression for the sidewash may easily be obtained by using the 
sidewash due to a horseshoe vortex in the same manner that the upwash due 
to a horseshoe vortex was used in obtaining equation ( 38 ). 


Approximation of bent lifting lines .- A bent lifting line can be approximated by a series 
of vortices of the form shown in figure 7 plus terms giving the effect of the bend in the center 
of the lifting line. Figure 8 shows a series of vortices of the form shown in figure 7 distrib- 
uted along a bent lifting line. The upwash from such a system of vortices is 
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(39) 
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where np is related to negative values of y and n2 is related to 
positive values of y. The first term on the right side represents the 
vortices for negative values of y; the second represents the vortices 
for positive values of y. The remaining terms take into account the 
effect of the bend in the center of the lifting line. Equation (39) is 
set up for a distribution of vortices as shown in figure 8. 

Equation (38) can be used to approximate a bent lifting line if X 
is replaced by Xp where Xp denotes x - xp. When this is done, the 
result is 


Equation (40) approximates a bent lifting line by adding the effects of a 
series of horseshoe vortices as shown in figure 9(a) when the circulation 
is symmetrical about the center of the lifting line, and when the points 
( x i,yi) are chosen symmetrically about the center of the lifting line. 

When the circulation is not symmetrical and/or the points (xp,yp) are 
not chosen symmetrically about the center of the lifting line, equation (40) 
approximates a bent lifting line by a series of vortices as shown in 
figure 9(b). The downwash from an infinite line vortex of constant 
strength and infinite slope is zero; therefore, equation (i»0) should give 
a reasonable approximation to a bent lifting line. 


In the theoretical development exact and approximate formulas have 
been developed for the potential and the upwash due to thin wings. An 
expression for the potential in space due to a thin wing was derived from 
an expression for the partial derivative of the potential with respect to 
the free-stream direction. The upwash was found by differentiation of this 
expression. Formulas are derived for the upwash in the wake and on the 
plan form. The results of the exact expression for thin wings were used 
to find the potential and the upwash due to lifting lines, and the exact 
ex P r sssions for lifting lines were used to find approximate expressions 
for lifting lines. 

In actual calculations better approximations to lifting surfaces 
may be obtained by using several lifting lines or their approximations 
distributed along the chord. Another approach would be to distribute 
finite vortices over the plan form. Either of the preceding methods 
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may be used to find a good approximation to the downwash from a lifting 
surface, provided that the apex of the hyperbola formed by the inter- 
section of the Mach cone with the plane of the wing does not lie on the 
plan form. For cases where the apex is on the plan form there is a 
finite contribution along the Mach cone. This contribution may be 
evaluated by the same method that was used to evaluate the double inte- 
gral in equation (15). 

It should be remembered that this paper does not discuss the 
effects of thickness; however, these effects may be evaluated rather 
simply by use of source distributions. 


The results of the theoretical development can be applied to a number 
of problems. A few of these problems will be considered here. 


APPLICATIONS 


Potential and Upwash at Infinity 


When x approaches infinity equation (6) becomes 



Since 


A0(x,y,z) = 


/ A 0xp dx 
J L.E. 


integrating the preceding expression with respect to xp gives 



Equation (Ll) is the same equation as the equation for the potential at 
infinity behind an airfoil in subsonic flow and is a well-known result 
noted in references 1 and 2. 
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The upwash at infinity is given by 



-)A0te 




! hl (Y 2 + z 2 ) 2 

Integrating equation (42) by parts and replacing A^ijrg by F gives 


(42) 
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1_ 

2n 


hi 


dT 

1 3yl 

Y 2 + z 2 
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( 43 ) 


For wings r(hp) and r(b2) are zero, the first term of equation (43) 
becomes zero, and the remaining term gives the upwash at infinity. The 
same result is given in reference 6. 


The Upwash for a Flat Unswept Wing of Infinite Span 
The A0 X on the plan form of an infinite unswept wing is 


Since 


c*A0x]_ 

dyi 


2 ccV 

’ P 


equation ( 17 ) reduces to 

0 z ( x >y>°) = _aV 


For points in the wake, equation (l4) gives 

0z( x >Y>°) = 0 

The preceding results are the well-known equations for the upwash in the 
z = 0 plane for the two-dimensional, flat, unswept wing. 


Camber and Twist Necessary to Give Specified 
Lift Distributions 

Equation (17) may be used to find the camber and twist necessary to 
give a specified lift distribution. The local angle of attack is given 

hy 

= _0_Z 
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and the pressure coefficient is given by 


2A0x]_ 


P = - 


(45) 


The local angle of attack for a given lift distribution, derived from 
equations ( 17 ), (44), and (4-5), is 


dP 


■R 


“ - ^ f ™ - hffr 471 


k J MC YX 


(46) 


As an example, the local angle of attack of a rectangular wing with 
constant load is now considered. 

When the pressure distribution is a constant, equation (46) becomes 

a = -T ' r* 4 — dXi (47) 

k knJ m YX 1 


For the rectangular wing in the region unaffected by the tip, 
equation (47) becomes 


• ■ -¥ 


(48) 


In the region affected by the tip, equation (47) becomes 

flO 


a = -£ £ - _L. 
k 4ny 




x - xp) 2 - P 2 y x 2 


x+0y 


(x - x x ) 


dx-i 


(49) 


where the coordinates are now located at the intersection of the leading 
edge and the tip. (See fig. 10.) After the integration is performed 
equation (49) becomes 


P / + \g2Zf + e tan-l 


Py 


4« \ 2 


7* 2 - P 2 y 2, 


(50) 


Equation (50) agrees with the results in reference 10. 
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Upwash Close to the Trailing Edge 

The value of the vertical perturbation velocity immediately behind 
the trailing edge can be obtained from equation (17)« 

In this application the concept of cancellation of pressure is 
utilized. It is assumed that the wing extends past the trailing edge. 
Since the pressure must be zero through the wake, this pressure behind 
the trailing edge must be canceled. The boundary condition^ in the wake 
are then satisfied by the wing pressure extended past the trailing edge 
and the cancellation pressure. The upwash is therefore made up of the 
effect of the pressure on the wing, its extension behind the wing, and 
the cancellation of the pressure behind the trailing edge. 

For trailing edges which are perpendicular to the free-stream direc- 
tion the application of equation (17) (using the cancellation-of -pressure 
concept) to the upwash immediately behind the trailing edge leads directly 
to the relation 


(0z)ij»g -o'te^ + 2^0x)fjrg 


(51) 


*here is the local angle of attack of the wing at the trailing edge. 

Che result given by equation (51) was obtained in references 1 and 2 from 
i physical consideration. When the trailing edge has a slope m, equation 
[ 17 ) yields 


(0z) 


TE 


~ 0 TE V 


(A0 x )te \/p 2m2 " 1 


2m 


(52) 


Equation ( 52 ) is also obtained in reference 1. 


Upwash in the Wake behind a Rectangular Wing 
Pitching about Its Leading Edge 

As an example to illustrate downwash calculations the upwash in the 
plane of the wing from a rectangular-pitching-wing is now calculated. 
(See fig. 10 for axes used in analysis.) For stationary axes a pitching 
wing moves in the arc of a circle; therefore, for axes fixed in the wing 
a pitching wing in linearized flow can be replaced by a wing which has a 
local angle of attack that varies linearly in the free-stream direction. 
It is assumed that the wake remains in the z = 0 plane and that the 
rolling-up of the trailing vortices can be neglected. 
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The wake behind a rectangular wing may be divided into two regions. 
One region is not affected by the tip or tips and is therefore a region 
which has a two-dimensional flow. The other region is the region affected 
by the tip or tips; in this region the flow is three-dimensional. 


There is no upwash in the wake behind a two-dimensional unswept, 
pitching wing; therefore, only the region affected by the tips of the 
wing has upwash. Since the analysis is based on linearized flow, the 
two tip effects add directly and therefore only one tip need be treated. 

The upwash in the z = 0 plane due to one tip of a rectangular 
pitching wing is shown in figures 10 and 11. The velocity potential and 
the pressure on a rectangular wing pitching about the leading edge were 
obtained by transforming the expressions given in reference 12 for the 
velocity potential and the pressure on a rectangular wing pitching about 
the half-chord line. 

The upwash at the trailing edge was found by using equation (51). 

The upwash close to the trailing edge was found by use of equation (17). 

The evaluations of the integrals were performed numerically. Equation (17) 
was derived to find the downwash on the plan form, but it may be used 
in the wake by applying the cancellation-of -pressure concept. 

A lifting line was used for calculating the upwash 2 or more chords 
behind the trailing edge inasmuch as this method was found to be accurate 
to two decimal places in this region. 

The upwash at infinity was found by using equation (^3). At infinity 
the upwash in the z = 0 plane for (3 = C = 1 is given by the following 
expressions: 

For y > 0 


0 Z = 




(53a) 


for -1 < y < 0 



(53b) 


for y < -1 


0Z - 



(53c) 
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Approximation of the Upwash from a 
Rectangular Pitching Wing 

The calculation of the downwash from the rectangular pitching wing 
together with results of reference 6 indicates that, for regions some 
distance behind the trailing edge, lifting lines give very good approxi- 
mations to the actual flow fields. In many cases it is simpler to use 
the approximation given by equation ( 38 ) than to use the exact lifting- 
line expression. The accuracy of equation ( 38 ) was investigated by cal- 
culating the upwash from an -unbent lifting line and comparing the results 
with the results obtained using equation ( 38 ). The unbent lifting line 
used was one that approximates the flow from the pitching rectangular 
wing. 


In figure 12 the values of the upwash from the exact and the approxi- 
mate lifting lines are plotted. Ten equally spaced points across the tip 
region were used in the approximation. It can be seen that for this case 
in most regions the approximation yields results which are as reliable as 
the exact values. If greater accuracy is required, it should be remem- 
bered that the accuracy of the approximation increases as the number of 
points is increased. 

Approximation to the Upwash from a 
Bent Lifting Line 

The approximation to the upwash from a bent lifting line given. by 
equation’ (39) was used to find the downwash along the line z = 0, y = 0 
behind a triangular wing as shown in figure 13 . Figure 13 shows the 
values given by equation (39) for 10 points across the semispan, the exact 
lifting line, and the exact linearized solution. The values for the exact 
lifting line for this case were taken from reference 6 and the values of 
the exact linearized solution were taken from reference 3* 

Figure 13 shows good agreement between the approximate and exact 
bent lifting lines. In general, the accuracy of the approximation will 
increase as the number of points used in the approximation is increased; 
therefore, it is expected that the difference between the exact and 
approximate expressions given in figure 13 could be decreased by increas- 
ing the number of points across the semispan. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., April 12, 1950 
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APPENDIX 

DIFFERENTIATION UNDER THE INTEGRAL SIGN 
OF THE EXPRESSION FOR 0 

The fact that equation (6) may be differentiated with respect to z 
under the integral sign without consideration of the variable limits can 
be proved as follows: 

If it is assumed that the right side of equation (6) can be 
differentiated by differentiating only the integrand, then 


: 2 + 


z 


If the preceding equation is correct, the potential may be expressed by 


*■-5 


f rr A0 X L( y2 _ z2 ) 


( Y 2 


3 2 z 2 (y 2 + z 2 )] 


X 


0 2 R3 


dx l d Yl 


(Al) 


, i f r z rr A ^x K y2 - *- 2 ) r2 + p 2 *- 2 fr 2 + x2 )] x , \ 

0 = i / / — 2J= dX! dy-L dX (A2) 

2 * (Y2 + x2) 2 [r(X)]3 


since 


0 = / 0z(x,y,X)dX 


Interchanging the order of integration of equation (A2) yields 
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The result of performing the first integration in equation (A3) is 


0 1 , f ^ dxj. 

2n JJ r (y2 + z 2)r 


d-yi 


(A4) 


Equation (A4) is the same as equation (6); therefore, it has been shown 
that equation (6) may be differentiated with respect to z under the 
integral sign without regard to the variable limits. Similar proof for 
the other velocity components can also be obtained. 
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Figure 1.- The intersection of forward Mach cone from point (x,y, z) with 

arbitrary plan form. 
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(b) Arbitrary point above wake. 
Figure 1.- Concluded. 
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Figure 2.- Regions of integration for a point affected by the wake. 
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Figure V.- Top view of regions of integration used in finding the down- 

wash on the plan form . 
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(a) Bent lifting line. 



(b) Vortex of constant strength. 
Figure 5 *- Lifting line of constant slope. 
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(a) Approximation that leads to equation (36). 



(b) Approximation that leads to equation (37). 



(c) Approximation that leads to equation (38). 
Figure 6.- Approximations of the 


span loading by series of rectangles. 
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Figure 7.- Finite vortex of the form used to approximate a "bent lifting 

line . 



Figure 8.- Bent lifting line approximated by a series of finite vortices. 
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(b) Unsymmetrical loading and/or unsymmetrically distributed 

points • 


Figure 9*- Bent lifting lines approximated by a series of horseshoe 

vortices . 
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Figure 10.- Contour plot of the downwash from a pitching rectangular 

wing in z = 0 plane. 



Figure 11.- Downwash behind a pitching rectangular wing in z = 0 plane. 
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Figure 12.- Downwash from a lifting line and its approximation by a series 
of horseshoe vortices at ten equally spaced points. (Lifting line 
approximates a rectangular pitching wing.) 



Figure 13.- Downwash determined by exact and approximate lifting lines 
alone the line z — 0, y = 0 behind a triangular wing with aspect 
ratio of 3.2. M = \p2.. 
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